Quadratic Dynamical Decoupling: Universality Proof and Error Analysis 
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We prove the universality of the generalized QDDjvjA^ (quadratic dynamical decoupling) pulse sequence 
for near-optimal suppression of general single-qubit decoherence. Earlier work showed numerically that this 
dynamical decoupling sequence, which consists of an inner Uhrig DD (UDD) and outer UDD sequence using 
JVi and Na pulses respectively, can eliminate decoherence to 0(T N ) using 0(N 2 ) unequally spaced "ideal" 
(zero-width) pulses, where T is the total evolution time and N — ATj = AT 2 . A proof of the universality of QDD 
has been given for even ATj . Here we give a general universality proof of QDD for arbitrary N\ and N2 . As 
in earlier proofs, our result holds for arbitrary bounded environments. Furthermore, we explore the single-axis 
(polarization) error suppression abilities of the inner and outer UDD sequences. We analyze both the single-axis 
QDD performance and how the overall performance of QDD depends on the single-axis errors. We identify 
various performance effects related to the parities and relative magnitudes of Ni and N2 . We prove that using 
QDD NlN2 decoherence can always be eliminated to o(T xaia{Ni ' N2} ). 
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I. INTRODUCTION 

The inevitable coupling between a quantum system and its 
environment, or bath, typically results in decoherence It 
is essential in quantum information processing (QIP) to find 
protection against decoherence, as it leads to computational 
errors which can quickly eliminate any quantum advantage 
IH [H. A powerful technique that can be used to this end, 
adapted from nuclear magnetic resonance (NMR) refocusing 
techniques developed since the discovery of the spin echo ef- 
fect UH], is dynamical decoupling (DD) @]. It mitigates the 
unwanted system-bath interactions through the application of 
a sequence of short and strong pulses, acting purely on the 
system. DD is an open-loop technique which works when the 
bath is non-Markovian and bypasses the need for mea- 
surement or feedback, in contrast to closed-loop quantum er- 
ror correction (QEC) 0]. However, while QEC can be made 
fault tolerant 10, H], it is unlikely that this holds for DD as 
a stand-alone method, or that it holds for any other purely 
open-loop method, for that matter fioll . This notwithstanding, 
DD can significantly improve the performance of fault toler- 
ant QEC when the two methods are combined ifTHl . 

DD was first introduced into QIP in order to pres erve 
single-qubit coherence within the spin-boson model Ill2l4l4ll . 
It was soon generalized via a dynamical group symmetriza- 
tion framework to preserving the states of open quantum sys- 
tems interacting with arbitrary (but bounded) environments 
H5I [Trill . These early DD schemes work to a given low or- 
der in time-dependent perturbation theory (e.g., the Magnus 
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or Dyson expansions lfl7ll ). Namely, the effective system-bath 
interaction following a DD pulse sequence lasting for a total 
time T only contains terms of order T N+1 and higher, where 
typically 7Y was 1 for the early DD schemes. For general N 
this is called iV th -order decoupling. Concatenated DD (CDD) 
Hill , where a given pulse sequence is recursively embedded 
into itself, was the first explicit scheme capable of achieving 
arbitrary order decoupling, i.e., CDD allows N to be tuned at 
will IU9f|. C DD has been amply tested in recent experimental 
studies ll20M24ll . and demonstrated to be fairly robust against 
pulse imperfections. However, the number of pulses CDD re- 
quires grows exponentially with N. In order to implement 
scalable QIP it is desirable to design efficient DD schemes 
which have as few pulses as possible. 

For the one qubit pure dephasing spin-boson model, Uhrig 
discovered a DD sequence (UDD) which is optimal in the 
sense that it achieves Ath order decoupling with the smallest 
possible number of pulses, N or N + 1, depending on whether 
N is even or odd j25ll . The key difference compared to other 
DD schemes is that in UDD the pulses are applied at non- 
uniform intervals. This optimal pulse sequence had also been 
noticed in lErjll for N < 5. A scheme to protect a known two- 
qubit entangled state using UDD was given in Ref . |27j] . UDD 
was conjectured to be model-independent ("universal") with 
an analytical verification up to N = 9 H and N = 14 
A general proof of universality of the UDD sequence was first 
given in IB0I1 (see also Ref. Bill for an alternative proof). The 
performance of the UDD sequence was the subj ect of a wide 
range of recent experimental studies Il2ll 132 - 13611 . An interest- 
ing application was to the enhancement of magnetic resonance 
imaging of structured materials such as tissue [37]. However, 
one conclusion from some of these studies is that the superior 
convergence of UDD compared to CDD comes at the expense 
of lack of robustness to pulse imperfections. It is possible 
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that recent theoretical pulse shaping developments |T3SL [39tl . 
designed to replace ideal, instantaneous pulses with realistic 
pulses of finite duration and amplitude while maintaining the 
suppression properties of UDD, will lead to improved experi- 
mental robustness. 

The UDD sequence is effective not only against pure de- 
phasing but also against longitudinal relaxation of a qubit cou- 
pled to an arbitrary bounded environment [f30h - That is, UDD 
efficiently suppresses pairs of single-axis errors. However, it 
cannot overcome general, three-axis qubit decoherence. The 
reason is that UDD uses a single pulse type (e.g., pulses along 
the x-axis of the qubit Bloch sphere), and system-bath inter- 
actions which commute with this pulse type are unaffected by 
the sequence. 

Combining orthogonal single-axis CDD and UDD se- 
quences (CUDD) reduces the number of control pulses re- 
quired for the suppression of general single-qubit decoher- 
ence compared to two-axis CDD lEol] , However, CUDD 
still requires an exponential number of pulses. This scaling 
problem was overcome with the introduction of the quadratic 
DD (QDD) sequence by West et ai, which nests two UDD 
sequences with different pulse types and different numbers 
of pulses Ni and N2 04111 . We denote this sequence by 
QDDA^jVa, where N% and N2 are the numbers of pulses of 
the inner and outer UDD sequences, respectively. QDDjv.jv 
(where the inner and outer UDD sequences have the same de- 
coupling order) was conjectured to suppress arbitrary qubit- 
bath coupling to order N by using 0(N 2 ) pulses, an expo- 
nential improvement over all prev iously known DD schemes 
for general qubit decoherence II4T[]. This conjecture was based 
on numerical studies for N < 6 14111 . and these were recently 
extended to N < 24 J42ll . in support of the conjecture. An 
early argument for the universality and performance of QDD 
(which below we refer to as "validity of QDD"), based on an 
extension of UDD to analytically time-dependent Hamiltoni- 
ans l43ll . fell short of a proof since the effective Hamiltonian 
resulting from the inner UDD sequences in QDD is not ana- 
lytic. 

The problem of finding a proof of the validity of QDD was 
first successfully addressed by Wang & Liu 14411 . though not 
in complete generality, as we explain below. In fact Ref. ll44tl 
considered the more general problem of protecting a set of 
qubits or multilevel systems against arbitrary system-bath in- 
teractions, using a nested UDD (NUDD) scheme, a gener- 
alization of QDD to multiple nested UDD sequences. This 
problem was also studied, for two qubits, by Mukhtar et 
ai, whose numerical results showed, for their specific choice 
of pulse operators, that the ordering of the nested UDD se- 
quences impacts performance l45ll . Wang & Liu's proof is 
based on the idea of using mutually orthogonal operation 
(MOOS) sets — mutually commuting or anti-commuting uni- 
tary Hermitian system operators — as control pulses 114411 (the 
ordering effect observed in Ref. l45ll disappears when using 
MOOS sets). As in QDD, the decoupling orders of the nested 
UDD sequences in NUDD can be different, so that different 
error types can be removed to different orders. Wang & Liu 
proved the validity of the general QDD/NUDD scheme when 
the order of all inner UDD sequences is even (the order of the 



outermost sequence can be even or odd) ll44Tl . Their proof is 
based on MOOS set preservation, and does not apply to QDD, 
or more generally NUDD, when the order of at least one of 
the inner UDD sequences is odd. In addition, Wang & Liu 
pointed out that there are QDDn 1 n 2 examples showing that 
the outer level UDD sequence does not work "as expected" 
(i.e., does not suppress errors to its order) if N\ is odd and 
N\ < N2- Thus their proof left the actual suppression order 
of QDD/NUDD with odd order UDD at the inner levels as an 
open question. 

This problem was addressed numerically in Ref. l42tl . 
which studied the performance of ODD^r^ for all three 
single-axis errors. The numerical results show that the sup- 
pression ability of the outer UDD sequence is indeed hindered 
by the inner UDD sequence if N\ is odd and, surprisingly, 
smaller than half of the order of the outer level UDD se- 
quence. Moreover, Ref. ll42tl reported that the suppression or- 
der of the system-bath interaction which anti-commutes with 
the pulses of both the inner and outer sequences depends on 
the parities of both N\ and A^. 

In this work we provide a complete proof of the validity of 
QDDjViJVa- 1° particular, we also prove the case of odd Ni 
left open in Ref. [|44h . Moreover, we analyze the single-axis 
error suppression abilities of both the inner and outer UDD 
sequences, and thus provide analytical bounds in support of 
the numerical results of Ref. ll42ll . 

We show that the single-axis error which anti-commutes 
with the pulses of the inner sequence but commutes with those 
of the outer sequence is always suppressed to the expected or- 
der (A^). The suppression of the two other single-axis errors 
(the one which commutes with the inner sequence pulses but 
anti-commutes with the outer sequence pulses, and the one 
which anti-commutes with both), is more subtle, and depends 
on the relative size and parity of N\ and N2 ■ 

Specifically, we show that when N± is even, ODD^^ al- 
ways achieves at least the expected decoupling order, irrespec- 
tive of the relative size of Ni and A^. However, when Ni is 
odd and Ni < N2 — 1, we show that the decoupling order 
of the error which commutes with the inner sequence pulses 
but anti-commutes with those of the outer sequence, is at least 
Ni + 1, smaller than the expected suppression order (A^). 
Nevertheless, for odd Ni and Ni > N 2 — 1, the outer UDD 
sequence always suppresses the error which commutes with 
the inner sequence pulses but anti-commutes with those of the 
outer sequence to the expected order (A^). 

One might expect that the error which anti-commutes with 
the pulses of both the inner and outer sequences can be sup- 
pressed by both sequences. In other words, one might expect 
this error to be removed at least up to order max[JVi, N2). 
However, we show that this expectation is fulfilled only when 
A^i is even. When Ni is odd, it determines the suppression 
order. However, interestingly, the parity of N2 also plays a 
role, namely, when it is odd the suppression order is one order 
higher than when N2 is even. 

Despite this complicated interplay between the orders of 
the inner and outer UDD sequences, resulting in the outer 
sequence not always achieving its expected decoupling order 
when A^i is odd, we show that, overall, QDDtv^ always sup- 
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presses all single-qubit errors at least to order min[iVx, N 2 ]. 

A complete summary of our results for the different single- 
axis suppression orders under QDD Xi n 2 is given in Table 
IIVI Our analytical results are in complete agreement with 
the numerical findings of Ref. J42j], but our proof method 
underestimates the suppression of of the error which com- 
mutes with the inner sequence pulses but anti-commutes with 
those of the outer sequence: for odd Ni we find a decou- 
pling order of min[iVx + 1,^2], while the numerical result 
is min[27Vi + l,iV 2 ] for N 1; N 2 < 24. Explaining this dis- 
crepancy is thus still an open problem. 

The structure of the paper is as follows. The model of 
general decoherence of one qubit in the presence of instan- 
taneous QDD pulses is defined in SecHJ The QDD theorem 
is stated there as well. We prove the QDD theorem in Sec. Hill 
and Sec.HV] A comparison between the numerical results of 
Ref. l42ll and our theoretical bounds is presented in Sec. [V] 
We conclude in Sec. [VI] The appendix provide additional 
technical details. 





inner oi 
even 


-der Ni 
odd 


outer 01 
even 


-der N 2 
odd 


no. of pulses 


Ni 


Ni + 1 


N 2 


Ni + 1 


no. of intervals 


Ni + 1 


Ni + 1 


N 2 + l 


N 2 + l 



TABLE I: Inner and outer sequence orders Ni and N2 vs the number 
of pulses and pulse intervals in the inner and outer sequences. 



N 2 = N 2 + 1 if N 2 odd. The additional pulse applied at 
the end of the sequence when N 2 is odd, is required in order 
to make the total number of X pulses-type even, so that the 
overall effect of the X-type pulses at the final time T will be 
to leave the qubit state unchanged. Note that it is the relative 
size of the pulse intervals that matters for error cancellations 
in UDD, not the precise pulse application times. Hence, the 
most relevant quantities for the outer level UDDjy 2 are the 
N 2 + 1 normalized UDDtv 2 pulse intervals (or the normalized 
QDDjVjjVa outer pulse intervals), 



II. SYSTEM-BATH MODEL AND THE QDD SEQUENCE 



A. General QDD 
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We model general decoherence on a single qubit via the 
following Hamiltonian: 

H = J Q I(g)B + J x a x ®B X + Jyo Y ®B Y + Jz°z ® B z , 

(1) 

where B\, A £ {0, X, Y, Z}, are arbitrary bath-operators with 
||i?^|| = 1 (the norm is the largest singular value), the Pauli 
matrices, a\, A 6 {X, Y, Z}, are the unwanted errors acting 
on the system qubit, and J\, A G {0, X, Y, Z}, are bounded 
coupling coefficients between the qubit and the bath. 

The QDDjViA^ pulse sequence is constructed by nesting a 
Z-type UDDjvj sequence, designed to eliminate the longi- 
tudinal relaxation errors ax <8> Bx and ay ® By up to or- 
der T Nl+1 by using Ni or N\ + 1 pulses, with an X-type 
UDD n 2 sequence, designed to eliminate the pure dephasing 
error a z <8> Bz up to order T N2+1 by using N 2 or N 2 + I 
pulses. The nesting order does not matter for our analysis, so 
without loss of generality we choose Z-type UDD Xl to be the 
inner sequence and X-type UDD x> 2 to be the outer sequence. 
We use the notation X and Z to denote control pulses, to dis- 
tinguish the same operators from the unwanted errors denoted 
by the Pauli matrices. We also sometimes use the notation 00 
for the 2x2 identity matrix /. 

The X-type UDDjy 2 pulses comprising the outer layer of 
the QDO^jVa sequence are applied at the original UDDjy 2 
timing with total evolution time T, tj = Tr/j where r/j is the 
normalized UDD timing (or the normalized QDDjViJVa outer 
sequence timing), 



rjj = sin 



J T 



(2) 



with j 



2{N 2 + 1) 

1,2,..., N 2 where N 2 = N 2 if N 2 even and 



where Tj = tj — tj-i is the actual pulse interval. 

The Z-type pulses of the inner level UDDa^ , applied from 
tj-i to tj, are executed at times 



3,k 



tj-i + Tj sin 2 



2(M + 1) 



with Ni + 1 pulse intervals 



Tj,k = tj.k ~ tj,k-i 



(4) 



(5) 



Even though adding an additional Z-type pulse to the end of 
each inner sequence with odd N\ is not required (since instead 
one can add just one additional Z-type pulse at the end of 
the QDDjViJVa sequence to ensure that the total number of Z 
pulses at the final time T is even), for simplicity of our later 
analysis, we let k = 1, 2, . . . Ni where Ni — N\ if Ni even 
and Ni = Ni + 1 if Ni odd. The corresponding normalized 
QDD jvi n 2 inner pulse timings are 



r)j,k 



T 



= Vj-i + Si sin 



I'TT 



2(Ni 



(6) 



with the normalized QDOAr^a inner pulse interval = 
Sj§k, where is the normalized UDD^ pulse interval and 
is the same function as Sj but with different decoupling order 
Ni . The first subindex stands for the outer interval while the 
second subindex stands for the inner interval. Moreover, by 
definition, we have 



T}j = Vj,Ni+l = Vj+1,0- 



(7) 



To summarize, the evolution operator at the final time T, at 
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the completion of the QDDtv^ sequence, is 



tions, 



-N 2 



with 



being the inner UDD n ± sequence evolution, and with U / be- 
ing the pulse-free evolution generated by H [Eq. ([TJ]. Table U 
summarizes how the number of pulses and pulse intervals in 
the inner and outer sequences depend on the inner and outer 
sequence orders N\ and N2. 
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B. Toggling frame 



Our QDD proof will be done in the toggling frame. Since 
our analysis is based on an expansion of powers of the total 
time T, most quantities we will deal with are functions of the 
normalized total time 1. 



The normalized control Hamiltonian with r\ = y is given 



by. 



N 2 N 2 + l Ni 

H c (v)^^[X^2 S (v-Vj) + Z £ (10) 

] = 1 j=l fe=l 



The normalized control evolution operator, 

C7 c (r?) =fexp[-i P H c ( V ')dr,'}, (11) 
Jo 

where T denotes time-ordering, is either / or Z in the odd j 
outer intervals, 



u c (v) = I, 
= z, 



[r)j,2£,'nj,2t+i) (12a) 

[%\2f+l,?7j,2£+2), (12b) 



while in the even j outer intervals, 



u c (v) = x, 

= Y, 



fe,2^,%',2«+i) (13a) 

%. 21+1 -.m. 21+2)- (13b) 



Accordingly, the normalized Hamiltonian in the toggling 
frame for the single-qubit general decoherence model, 

H(rj) = U e (itfHUM (14a) 
= foJol ® B Q + f x (T])J x a x ® B x (14b) 
+fy{ r i)Jy (J Y ®B Y + fz(v)Jz°z <8 B z , 



= fx(r])fz(v) 



[0,1,), (15a) 
[Vj-uVj), (l 5b ) 
[Vj,k-i,Vj,k), (15c) 
[Vj,k-i,Vj,k), (15d) 
(15e) 



unlike the single-qubit pure dephasing case, which has only 
two UDD modulation functions. Because the Z-type pulses 
on the inner levels anti-commute with the errors ax and ay 
and commute with az, the modulation functions f x {r\) and 
fy(rf) switch sign with the inner interval index k while f z {v) 
is constant inside each outer interval. On the other hand, the 
outer X-type pulses anti-commute with the errors a z and ay 
and commute with the error ax, so both f z (i]) and f y (ij) 
switch sign with the outer interval index j, while f x (77) doesn't 
depend on the outer index j. 

Each QDDjv^ modulation function f\{rj) can be sepa- 
rated naturally as 



fx(v) = f*(v)h(v) 



(16) 



where /a (77) describes the behaviour of f\ (77) inside each 
outer interval and fp(r]) describes the behaviour of fx (77) 
when the outer interval index j changes. In fact /s(»7) is 
identified as the normalized UDD n 2 modulation function and 
fa (r?) covers N% + 1 cycles of UDD ^ modulation functions 
with different durations. However, up to a scale factor f a (v) 
is the same function in each of these cycles. Therefore, instead 
of f a {v)> we use one cycle of the normalized UDD^ modu- 
lation function denoted as /„ (77) to denote the effective inner 
function of f\{rj). Likewise, since fs(r]) is constant inside 
any jth outer interval Sj, it can be viewed as a function of the 
outer interval j, and we replace fp(rf) by the notation fp(j). 
In particular, fp =z (j) = (— l)- 7-1 . Table [TT] lists the effective 
inner functions f a and outer functions fp for all QDDjviJv 2 
modulation functions fx and will be used in Sec. [Ill] 



fx 


(fa, f(j) 


fx 


(fx,fo) 


fv 
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TABLE II: The effective inner functions f a and outer functions fp of 
the normalized QDD jv 2 jv 2 modulation functions fx . Functions in the 
first column are the normalized QDDjvj n 2 modulation function and 
those in the second column are the normalized UDDjvj and UDDjv 2 
modulation functions respectively. 



has four different normalized ODD^^ modulation func- 



In the toggling frame, the unitary evolution operator which 
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Components \nth order QDDjViA^ coefficients: 




\i=X 


a\ n .. 


Ai=Y 


Un- 


Ai = Z 


Total # of ax and f x (rj) 


odd 


even 


even 


odd 


even 


odd 


( 1 ) Total # of cry and f y (n) 


even 


odd 


odd 


even 


even 


odd 


Total # of az and / z (77) 


even 


odd 


even 


odd 


odd 


even 


(2) Total # of effective inner integrand f x 


odd 


odd 


odd 


odd 


even 


even 


(3) Total # of effective outer integrand f z 


even 


even 


odd 


odd 


odd 


odd 



TABLE III: Number combinations of Pauli matrices (or modulation functions) for each error type (or QDDjvj jv 2 coefficients). For example, 
when An • • • Ai = X, there are two possibilities, represented in the two corresponding columns in the rows numbered (1): either there is an 
odd number of ax (and f x (rj)) along with an even number of both ay (and f y (rf)) and az (and f z (r])) , or there is an even number of ax (and 
f x (ri)) along with an odd number of both ay (and f y (rj)) and az (and fzijl)) . Consulting Table HH in the first case there is an odd number of 
inner integrand f x functions from f x {v) an d an even number of f x from f y {ri), so that the total number of f x is odd, as indicated in the first 
entry in row (2). Likewise, in the second case there is an odd number of outer integrand f z functions from f y (rj) and an odd number of f z 
from fz(rj), so that the total number of f z is even, as indicated in the second entry in row (3). 



contains a whole QDDjvj at 2 sequence at the final time T reads 



C. Error terms 



U(T) = TexpH / H(t)dt 



(17a) 



= Texp[-iT / H(r])dr]}. (17b) 
Jo 

We expand the evolution operator U(T) into the Dyson series 
of standard time dependent perturbation theory, 

oo 

U{T) = ^2(-iT) n J^a^ ® B^a Xn ... M , (18) 
"=° a„ 

where we use the shorthand notation 

E = E E - E • c«> 

X„ \ n £{0,X,Y,Z} A„_iG{0,X,y.Z} Aie{0,X,F,Z} 



and 



^ n ^ =n^' B i n) = n^ A =- (20) 



Finally, 



t = / dr,^... / d v ^ l[h e (v {e} ) (2D 

Jo Jo e=1 



a A„-A 



is a dimensionless constant we call the nth order normalized 
QDDat^ coefficient. These coefficients play a key role in 
the theory as it is their vanishing which dictates the decoupling 
properties of the QDD sequence. 

The subscript of a\ n represents a product of Pauli ma- 
trices, and we shall write A„ • • • Ai = A with A representing 
the result of the multiplication up to ±1, ±i. From Eqs. ( fT8l - 
(f2Tb . this subscript indicates not only its associated operator 
term a\ n . . . <j\ 1 <X> B\ n ■ ■ ■ B\ 1 but also its n ordered inte- 
grands, fx n ■ ■ ■ fx ± . Moreover, from Table HT1 one can also 
deduce the ordered set of effective inner and outer integrands 
for a given subscript of a\ rl ...\ 1 . 



(n) 

Every product of system operators, a\ = a\ n ■ ■ ■ a\ 1 can 
be either /, ax, &y or &z- Th e summands in the expansion 
( fl~8b of U can accordingly all be classified as belonging to one 

of four groups . If G {crx, cry^z} men tne correspond- 
ing summand in Eq. ( fT~8T > decoheres the system qubit. 

Definition 1. A single-axis error of order n and type A is 
the sum of all terms in Eq. ( 1181 ) with fixed and A G 
{X,Y,Z}. 

In the UDD case Eq. ( fT8l would include just one type of 
single-axis error [f30h - 

Due to the Pauli matrix identities aiOj = ie^kCk and 
of = /, which of the three possible errors a given product 

(n) 

o~\ becomes is uniquely determined by the parity of the total 
number of times each Pauli matrix appears in the product. In 
this sense there are only two possible ways in which each type 
of error can be generated. Take the error ax as an example. 
One way to generate ax is to have an odd number of ax op- 
erators which generates ax itself, along with an even number 
of ay, an even number of az, and arbitrary number of /. The 
other possibility is an odd number of ay with an odd number 
of az to generate ax, along with an even number of ax and 
arbitrary number of /. 

Note that for a given error a\, the parity of the total num- 
ber of times each modulation function appears in aA„..-Ai=A' s 
integrands f\ n ■ ■ ■ f\ 1 is also determined accordingly. For ex- 
ample, consider A n • • • Ai = Z. This can be the result of there 
being an even number of az [and f z (v)] along with an odd 
number of both ax [and f x (v)] and ay [and f y {rj)], a situa- 
tion summarized in the last column of the block numbered ( 1 ) 
in Table [Til] In this case, given Table [TTJ the total number of 
effective inner integrand functions f x contributed by f x [rf) is 
odd, as is the contribution of effective inner integrand func- 
tions f x from f y {rj), so the total number of effective inner in- 
tegrand functions f x is even. This situation is summarized by 
the last "even" entry in row (2) of TablelTTTl This table gives all 
possible parities of Pauli matrices (or modulation functions) 
for each type of error (or its associated nth order QDDjv^ 
coefficient). The parity of the total number of identity matri- 
ces I (modulation function /o) is irrelevant for the proof, so 
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is omitted from Table [TTT] With a given number combination 
of Pauli matrices (or modulation functions) and Table ITT1 one 
can determine the parity of the total number of effective inner 
and outer integrands as presented in rows (2) and (3) of Table 
Hill Table Hill will be referred to often during the proof. 

If all of the first A^th order QDD a\ coefficients vanish 
for a given A, namely a\ n ...\ 1= \ = for n < N, we say 
that the ODD^y^ scheme eliminates the error a\ to order 
N, i.e., the error a\ is 0(T N+1 ). Naively, one might expect 
the inner Z-type UDD n 1 sequence to eliminate both ax and 
ay errors to order N±, and the outer X-type UDD^ sequence 
to eliminate az errors and any remaining ay errors to order 
N2. The situation is in fact more subtle, and is summarized in 
the following QDD Theorem whose proof is provided in Sec. 
|III]and Sec.|IV] 

QDD Theorem 1. Assume that a single qubit is subject to 
the general decoherence model Eq. ([T). Then, under the 
QDDn 1 n 2 sequence Eq. ([8), all three types of single-axis 
errors of order n are guaranteed to be eliminated if n < 
min[iVi, N 2 ]. Higher order single axis errors are also elim- 
inated depending on the parities and relative magnitudes of 
Ni and N 2 , as detailed in Table \Iy\ the results of which re- 
main valid under any permutation of the labels X, Y, Z along 
with a corresponding label permutation in Eq. ([8]). 



we find that the error ay is suppressed to order max[A r i , N 2 ] 
if N 2 is even, and to order maxf^ + 1,-/V 2 ] if N 2 is odd. 
These results are all summarized in Table II VI 



III. SUPPRESSION OF LONGITUDINAL RELAXATION 

ax AND cry 

A general proof of the error suppression properties of UDD 
was first given by Yang & Liu, including for the suppression 
of longitudinal relaxation errors ax and ay lf3Qh . Wang & Liu 
first proved that the outer sequence does not interfere with the 
suppression abilities of the inner sequence with the DD pulses 
chosen as a MOOS set l44ll . In this section, we give an alter- 
native non-interference proof which shows explicitly that it is 
the inner Z-type UDD^ sequence that makes all longitudi- 
nal relaxation related QDD^jVa coefficients a\ n ...\ 1=ax ^ Y 
with n < Ni vanish, regardless of the details of the outer 
X-type UDDat 2 sequence. Moreover, we also show that the 
outer X-type UDDat 2 sequence, when the outer order N 2 is 
odd, eliminates the ay error to one additional order, i.e., to 
order JVi + 1. For precise details refer to Table llVI 



Single-axis 


Inner order 


Outer order 


Decoupling order 


error type 


N x 


^2 






arbitrary 


arbitrary 


Ni 




even 


even 


max[Ni,N-z] 


ay 


even 
odd 


odd 
even 


max[iVi +1,N 2 ] 
Ni 




odd 


odd 


Ni + 1 




even 


arbitrary 


N 2 


odd 


arbitrary 


min[7Vi + 1,JV 2 ] 



TABLE IV: Summary of single-axis error suppression. For each 
error type ax, the nth order QDD n 1 n 2 coefficients [Eq. J21H 
a\ n ---\ 1 =\ = Vn < N, where V is the decoupling order given 
in the last column. 



An immediate corollary of this Theorem is that the overall 
error suppression order of QDDat^ is min[iVi, 7V2]. This 
will be reflected in distance or fidelity measures for QDD, 
such as computed for UDD in Ref . 13111 . 

We shall prove Theorem Q] in two steps. First, in Sec. Hill 
we shall prove that for arbitrary values of Ni and N 2 the 
QDDjVi^ sequence eliminates the first Ni orders of ax and 
ay errors. Secondly, we shall prove that if N 2 is odd, an ad- 
ditional order of the ay error is eliminated, i.e., Ni + 1. We 
will not show any suppression of the az error in Sec. [TlTl 

Then, in Sec.|lV]we shall complete the analysis of the effect 
of the outer sequence, and show that the az error is suppressed 
to order N 2 if Ni is even. If Ni is odd, az is suppressed to 
order N 2 if Ni > N 2 - 1, and to order Ni + 1 if iVi < 
N 2 — l. Additionally, we show that if Ni is even, the ay error 
is suppressed to order A^2, which may be higher than the result 
of Sec. Iim alone. Combining the results of the two sections, 



A. The outer interval decomposition of a\ n ...\ 1 

We expect the inner Z-type UDDjVi sequences of 
QDDjVjjVa to suppress the errors ax and ay. Therefore, our 
strategy for evaluating a\ n ...\ 1 [Eq. ( 1211) 1 is to split each of 
its integrals into a sum of sub-integrals over the normalized 
outer intervals sj in Eq. (l3al ). In this way, each resulting seg- 
ment of a\ n ...\ 1 can be decomposed naturally into an inner 
part (which contains the action of the inner Z-type UDD^rJ 
times an outer part (which contains the action of the outer 
X-type UDDjy 2 sequence). The manner by which the inner 
Z-type UDDjvj sequences suppress longitudinal relaxations 
can then be easily extracted. 

As we show in Appendix fAj after this decomposition 
dA„ -Ai can be expressed as 

flAn-Al = * ln ( r n/a„ r n-l ■ • ■ /a 2 r l/ai) X 

{r/=0,*}J- 1 I 

t> out (r n f Pn r n ^ 1 ...fp 2 r 1 f l3l ). (22) 

with r n = *. This is just a compact way of writing multiple 
nested integrals and multiple summations, with a notation we 
explain next. 

First, f at and fp t are the effective inner and outer functions 
respectively of a\ n ...\ 1 's £th integrand f\ e . From TableHTl the 
effective inner (outer) function of the normalized QDDjviJV 2 
modulation functions will be either f x (f z ) or /o = 1, the 
normalized UDDjVi (UDDat 2 ) modulations functions in the 
generic ax (o'z) pure bit flip (dephasing) model. 

Second, the "inner output function" $ m generates all the 
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segments' inner parts via the following mapping, 



nf ae 

For example, 



* m (*/ Q2 0/ Ql )= / fa 2 d V 

Jo Jo 



1 r rf- 
(2) 



faidr) 



(i) 



integrals. More specifically, if the sub-integrals over any two 
adjacent variables rf^' and rf i+1 ^ are already located in time- 
ordered, different outer intervals, then the sub-integral over 
(23) 7 j{() j s not nes i; ec i inside the sub-integral over rf i+1 \ and its 
integration domain is the entire outer interval. In contrast, 
if the sub-integrals over any two adjacent variables ??W and 
7^+!) ^e in the same outer interval, it follows that their sub- 
integrals are nested due to time-ordering. 

(24) 



a term which appears in the expansion of a\ 2 \ 1 . 

From Eq. ( 1231 , one can see that ri determines how the in- 
tegral of 77 W relates to the integral of its adjacent variable 
For the inner part, the relationship between the in- 
tegrals of two adjacent variables i^ e+1 '> and rf^ is either in- 
dependent (r^ = *; they appear in separate integrals) or nested 
(ri = 0; they appear together in a time-ordered pair of inte- 
grals). 

Third, the outer output function <I> out generates all the seg- 
ments' outer parts via the following mapping, 



"Tlfol 



where is the j^th normalized outer interval for variable 
77W, and m indicates that n is the mth * in {r n r n -i . . . 7*1}, 
counting from r±. For r n fp n with r n = *, the upper limit 
_ j(£+i) _ 1 in Eq. d25l should be replaced by j( n *> = 
N2 + 1. For example, 

* TOt (*/fc * /fc0/fc) = E /ft(i (3) )^(3) x 

j(3) =2 

E /ftC7' (2) )»jw/ft(j (a) )*j(»). (26) 

a term which appears in the expansion of a\ 3 \ 2 \ 1 . 

From Eq. d25l ). r£ indicates the relationship between the 
outer intervals of two adjacent variables r/^ +1 ) and rf '. They 
can either be time-ordered, namely, in different outer intervals 
(ri = *), or in the same interval (ri = 0). 



Finally, 



E 



E E 



E 



(27) 



{m=0,*}?=i r n _ie{0,*}r„-2e{0,*} rie{0,*} 



includes all possible integration configurations for $ m 
and all possible summation configurations for $ out . 
Each such configuration is determined by a given set 
{*,r n -i,r n -2, ■ ■ -,ri}. 

Note that the integration pattern of the inner part determines 
the summation pattern of the outer part and vice versa. The re- 
lation between the inner part and its corresponding outer part 
comes from the time-ordering condition, j/W > t?^™ -1 ) > 
. ..Ty?) > T)^\ because a^„..-Ai comprises n time-ordered 



B. The inner parts $ ln and the outer parts $ out ofa\ n ...\ 1 

Consider the argument r„/ Aln r n _i/ Al „_ 1 r„_ 2 • ■ • 
of i> m or $ out , where fx can be a or j3 in Eq. ( |22] |. Define 
a cluster of /'s as a contiguous set of /'s connected only 
by 0. Different clusters are separated by *, For example, 

(Ub ) * (/m4 0/(U 3 ) * (/^2 0/mi )• where the parentheses indicate 
clusters. In this manner, each integration or summation con- 
figuration {*r„_ 1 r„_ 2 . . • ri} corresponds to a way in which 
a set of n functions is separated into clusters. 

Suppose that for a given configuration {*r„_ir„_2 ■ • • fi}, 
the mth inner cluster, counting from right to left, is 
/ap0/a p _i0 ■ • • ®fa q , which has p — q + 1 elements. Likewise, 
we have the mth outer cluster fp p ®fp p _ 1 ® . . . ®fp q - Applying 
the rule of Eq. ( 1231 to the mth inner cluster, or the rule of 
Eq. d25l ) to the ?nth outer cluster, we then have, respectively 



/o,0-.-0/a,(*) 



drj^ I dr) 
Jo 

■*a„a p _i...Q„ 



(P-1) 



(28a) 

l=q 

(28b) 



jm+l-l 



j m —va t—q 



where if p = n, namely the mth group is the last group (count- 
ing from right to left), the upper limit j m +i — 1 should be 
replaced by j m = N2 + 1. Also note that in Eq. ( I28cb we 
have replaced [according to the notation of Eq. d25b l by 
the cluster index j m . 

Now recall that the outer effective function fp e (j) is either 
f = lor/ 2 (j) = (-ly- 1 . Therefore, if TIL, /A (J'™) 
contains an odd number of f z (j m ), we have fpt (jm) = 
(_l)im-i ) otherwise ]J p e=q ffs e (j m ) = 1. 

Note that the nested integral u ap a p -i...a q m Eq. ( I28ab is 
just the (p ~ q + l)th order normalized UDD^ coefficient 
for the generic ax pure bit flip model, because the effective 
integrands f at , either f x or /o, are the normalized LTDD^ 
modulations functions. 

We have now assembled the tools to perform the summation 
implied in Eq. d22l . which is the result of the outer interval 
decomposition. Different clusters, each of which is given in 
Eqs. d28al > or ( |28cl i, are simply multiplied. To illustrate this, 
the second order normalized QDDyv! n 2 coefficients a\ 1 \ 2 are 
listed in Table[V] 
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Error Type 


2nd order normalized QDD n ± n 2 coefficients 


ox 


a x0 = U x U Ei=2 Sj Ep=l S P + U xQ Ej=l s 3 
a 0x = U U x S 3 Ep=l S P + U 0x Ej-1 S j 

a zy = uo Ux e^c-i)^ 1 ^ Ei=i(-iy-% + «0« Egi +1 4 


Oy 


a y o = ^uoE J = 2 (-!) J s i Ep=i s P + u x0 Ej=i (-1) J «j 
a 0y = u x E^ 2 +1 *j Eg("irS + "o, Ef=i +1 (-l) J '-S- 
a M = u u x Zpt'i-iy^Sj Eg s P + u 0x E^rC-lP'- 1 ^ 

^ = E^ 2 2 +1 *s Z P r,S(-i) p -\ + «*o Efirc-ip'" 1 ^ 




Y^ v 2 + 1/ iyi'-l V^P<J i v^a+l/ ivi-1 2 

a z0 =woUoEj=2 l^ 1 ) s i Ep=i s p + u oo Ej=i s j 

= «o uo eS 2 + S Eg(-ir\ + «oo Eli'Hr 1 ^ 2 

O-xy — ^x ^x 2 ^^p=l( *^)^ "'sa; = l ( ^/ 
— ^x ^x j — ? ( ^a;^ /* > 7 — 1 ( *^)^ 



TABLE V: The outer decomposition form of oajAj . For n = 2 we have r2ri = {**} or rir\ = {*0}. The first summand in each line is the 
result of the {**} expansion, the second is the result of the {*0} expansion. 



The following lemmas relate the normalized QDD and 
UDD coefficients. They are easily concluded from Eq. ( 128a| ). 

Consider a configuration {*r„_ 1 r„_ 2 • ■ • r i} with m *'s. 
Correspondingly there are m clusters. Each cluster has associ- 
ated with it a normalized UDD n 1 coefficient of order n' equal 
to the number of elements (/'s) in the cluster, and 1 < n' < n. 
The sum of all the orders is n. Thus: 

Lemma 1. Consider a configuration {*r n _ir n _2 • • • with 
m * 's. The corresponding inner part $ ln ofa\ n [Eq. (1221 ) ] 
is composed of m normalized UDDx 1 coefficients whose in- 
tegrands are the effective inner ones of a\ n . ..x^ and the sum 
of whose orders is equal to n. 

In addition, all of the first nth order UDD n 1 coefficients, 
but not order n + 1 and above, appear in any given nth order 
QDDjVjAT;, coefficient, i.e., 

Lemma 2. The only UDD n ± coefficients which can appear in 
all the inner parts <I> m of the nth order QDDn 1 n 2 coefficient 
a\ n ...\ 1 are those whose orders are between 1 and n. 

C. The first Ni vanishing orders of the single-axis ax and cry 
error due to the inner Z-type UDD n x sequences 

From the second row of Table [Till one can see that the to- 
tal number of / x 's in the effective inner integrands of the co- 
efficients a\ n ...\ 1= x and a\ n ...\ 1= y is odd. Accordingly, no 
matter how one divides the inner integrands into clusters, there 
will always be at least one cluster which has an odd number of 
f x . Then, from LemmaQ] it follows that all the inner parts $ m 
°f Ox„...Ai=X an d a X n —\i=y contain one or more UDDjVi 
coefficients with an odd number of f x in the integrands. Re- 
call that UDD^ coefficients ux m ...Ai=x> i- e -> those associ- 
ated with the error ax , contain an odd number of f x in their 
integrands. Therefore, we have 

Lemma 3. After outer interval decomposition, all the in- 
ner parts of the nth order QDDn x n 2 coefficients a\ n ...\ 1 —x 



and a\ n ...x 1 —Y contain one or more UDDn 1 coefficients 
u\ m ...\ 1= x, where m < n. 

Now recall: 

Lemma 4. (Yang & Liu [3(\l ) The UDDn 1 coefficients 
u\ m ...\i=x = when m < N\. 

It follows from the last two lemmas that all QDDa^a^ co- 
efficients associated with longitudinal relaxation a\ n ...\ 1= x 
or a\ n ...\ 1= y with n < Ni vanish. Physically, it is clearly the 
inner Z-type UDD sequence that is responsible for elimi- 
nating the single-axis errors ax §S Bx and ay ® By up to 
order T Nl+1 . The effect of the outer X-type UDD^ se- 
quence is entirely contained in the outer output function <E> out 
in Eq. d22l . and consequently does not interfere with the elim- 
ination ability of the inner level control Z-type UDD^, in 
agreement with l44ll . 

Fromrow 2 of TablelllTl unlike a\ n ■■■\ 1 =x.y, all a\ n ...\ 1= z 
contain an even number of effective inner functions f x . Ac- 
cordingly, Lemma [3] does not apply to a\ n ...\ 1= z and there- 
fore, the argument that all the inner output functions <I> m of 
a\ n ...\ 1 are removed by the Z-type UDD^ sequences cannot 
be applied to the pure dephasing terms. This is, of course, due 
to the fact that the inner Z-type sequence commutes with the 
the pure dephasing error az C8) Bz- Instead, this error will be 
suppressed by the outer X-type UDDtv 2 sequence. 

Note that the outer output functions <I> out of a\ n ...\ 1 
[Eq. (1221)1. which contain the effect of the outer X-type 
UDDat 2 sequence, are expressed in terms of multiple time- 
ordered summations [Eq. ( 128cH . which are not easily ana- 
lyzed using the current method. Therefore, in order to demon- 
strate the suppression of the pure dephasing error az <8> Bz, 
in Sec.|IV]we shall deal directly with a\ n .\ 1= z, rather than 
a separation into inner and outer parts as we have done in this 
section. 
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D. One more order of suppression for the single-axis cry error 
due to the outer X-type UDDjv 2 sequence when N2 is odd 

In the previous subsection we proved that a\ n ...\ 1= Y = 
when n < Ni. Now we shall show that also a\ N +1 ...Xt=Y 
vanishes, due to the outer level X-type UDD jy 2 sequence, for 
odd N2- Essentially, as we now explain in detail, this sequence 
is responsible for eliminating one remaining term in the ex- 
pansion of a AjVl+1 ...Ai=y- 

According to Lemma|2l as applied to a\ N +1 ...\ 1= y, the 
only UDD jvi coefficients which can appear are those with or- 
der at most Ni + 1 . According to Lemma [4] the first N% of 
these UDD coefficients vanish. The only UDD coefficient (in 
o-\ Ni+1 ...\ 1 =y) regarding which at this point we have no infor- 
mation is the Ni + 1th, and indeed, it may be nonvanishing. 
Using the mapping Eq. ( I28cl ). we therefore have 



"A 



JVi + l 



...Ai=y 



N 2 + l iVi + 1 

3 = 1 t=l 



N 1 - 



(29) 

We now show that this vanishes due to the outer part. 

According to the third row of Table [Till all a\ n ...\ 1= Y 
contain an odd number of effective outer functions f z (j) = 

(-1) 7 ' -1 . Consequently, we have T]fci +1 /ft 0') = 
which simplifies the outer part in Eq. d29l to 



N 2 +i 

£(- 



-1 JVi+1 
b 3 



(30) 



Note that the UDD pulse intervals are time-symmetric (for the 
proof see AppendixlBt. Therefore, the UDD^ outer intervals 



satisfy 



SN 2 +2-j 



(31) 



If the decoupling order of the outer UDD sequence N2 is odd 
then j and N2 + 2 — j have opposite parities. Accordingly, 



( _ 1)3 _i = (-iy-i( SN2+2 _^^ (32) 



/ ■\\N 2 +2-j-l/„ \Ni+l 



Thus, when N2 is odd the outer part Eq. (f30b vanishes due 
to the mutual cancellation of terms with equal magnitude but 
opposite sign. 

This concludes our proof of the error suppression of ax and 
<jy errors to order JVi, and of ay to order N\ + 1 when X2 
is odd. This confirms row one of Table II V I and row two of the 
same Table, disregarding for now N2 in the last column. In 
the next Section we set out to complete the proof and confirm 
all claims made in Table II VI 



IV. SUPPRESSION OF THE PURE DEPHASING ERROR 

In this section we focus on the suppression of the pure de- 
phasing error az by the outer X-type sequence, and also show 



that <jy can be additionally suppressed by the outer sequence 
to order N2 when N\ is even. 

To do so, we shall show that if JVi is even the inner Z- 
type UDD jv x sequence does not hinder the suppression abil- 
ity of the Y and Z-type errors by the outer A-type UDD jy 2 
sequence. For odd Ai we cannot conclude that the inner se- 
quence does not hinder the outer sequence. However, if N\ is 
odd, our method does show that the outer sequence suppresses 
oz at least to order min[Ai + 1, A2]. 



A. Linear change of variables 

To avoid having to analytically integrate a multiple nested 
integral with step functions as integrands such as ax n ...\ 1 , we 
adapt the approach of Refs. ||30[ 14411 . which avoids integrating 
step functions directly but still manages to show a\ n ...\ 1= z = 
up to a certain order. 

First, we make an appropriate variable transformation from 
rj G [0, 1) to 9 G [0, 7r), with the result that the outer pulse 
intervals are all equal. This is required to make the modulation 
functions f x , f y , f z , and /o (possible integrands that can occur 
in a\ n ...\ 1 ) become periodic functions so that each of their 
Fourier expansions is either a Fourier sine or Fourier cosine 
series. 

The variable transformation introduced by |44| to tackle 
the QDDatjJVj, sequence is to apply the corresponding linear 
transformation to each outer pulse interval [r/j-i , r]j) with du- 
ration Si 



N, 



V - Vj-i 



(3 - j> 
A 2 + l ' 



The timing of the outer A-type pulses becomes 



9 = J7T 
3 N2 + I 



(33) 



(34) 



so that f z (9) becomes a periodic function with period of N 2 2 ^ 1 , 

f z (9) = i-iy- 1 [Oj-uOj) 

00 

= J2 d k sin[(2A; + l)(A 2 + 1)0], (35) 

fc=0 

where the second equality is the Fourier sine-series expansion, 
and dt 



l k ~ (2A:+l)7r p 

When we apply the piecewise linear transformation d33b to 
the inner pulse timings r\j^ [Eq. ©] the UDD^ structure is 
preserved 



N 2 + l V2(Ax + l) 



k'TT 



(36) 



In fact all the inner pulse sequences become identical as they 
have the same total duration N * +i . It follows that f x (0) = 
(— l) fc_1 within [6j k-i, 0j,k) is a periodic function with pe- 
riod of tv^t-t. 
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The parity of the decoupling order Ni of the inner UDD^V! 
sequence determines whether f x (9) is even or odd inside each 
outer interval (Appendix 0. Inside each outer interval the 
parity of f x (0) equals that of Ni. Hence, we have 



fx(0) 



E d% cos[2fc(A 2 + 1)6»] Ni even 



fc=0 



(37) 



E d x k sin[2fc(A 2 + 1)0] Nt odd 



U-=i 



The relation f y (0) = f z {9)f x {8), Eqs. (03) and (07), and 
the product-to-sum rules of the trigonometric functions, 



sin a sin b 



[cos (a — b) — cos (a + b)] , (38a) 



cos a sin b = — [sin (a + b) — sin (a — b)] , (38b) 
cos a cos b = - [cos (a + 6) + cos (a — 6)] , (38c) 
yield the following Fourier expansions of f y (9) 

' oo 

E d\ sin[(2fc + l)(iV 2 + 1)0] iVi even 



fy(0) 



k=0 



(39) 



E^ cos[(2fc + 1)(7Y 2 + 1)6>] iVi odd 



fc=0 



Note that while the Fourier expansion coefficients in the 
even and odd cases are in fact different, we use the notation 
df or d v k for both since the exact values of these coefficients 
are irrelevant for our proof. 

It follows from Eq. (03) that d V = iY ^ ti s J dO = G(0)d6, 
where G(0) is the step function whose step heights are pro- 
portional to the ODDITY;, outer intervals, 



G(0) 



N 2 + l 



0e [0,-1,0,) (40a) 



E E 9k, q sm[^k)(N 2 + l)d + qd],(40b) 

k=0 q=—l,l 



as shown in AppendixlDl 

With Eqs. (05), d37)-(l40bl. and f = 1, the nth order 
QDDA^TVa coefficients (01) can be rewritten as 



d0n / dfl„_l 

o Ja 



'Milium 



with / A ,(0 f ) = G{6i)f Xl {e t ), where 



/o = E Yl 9k, q sm[2k(N 2 + l)9- 

k=0 s=-l,l 



(41) 



(42) 



/« = X) E 4 g cos[(2fc + l)(7y 2 + l)0 + g0].(43) 
fe=o «=-l,l 



When the inner decoupling order Ni is even, 

oo 

h = E E dl q sm[2k(N 2 + l)e + q ei (44) 

fc^o g=-l,l 

oo 

fv = E E < g cos[(2fc+l)(iV 2 + 1)0 + ^0] ,(45) 

fc=0 g=— 1,1 

while if it is odd, 

oo 

/» = E E d% iq cos[2k(N 2 + l)9 + q9], (46) 

k=0 9=-l,l 
oo 

= E E < g sin[(2fc + l)(iV 2 + l)0 + g0].(47) 

fc=0 g=— 1,1 

Observe that all the integrands of aA n ---Ai are composed of 
sums of either purely cosine functions or purely sine func- 
tions, i.e., none of the integrands is a mixed sum. This fact is 
key to our ability to perform the nested integral, as we show 
next. 



B. Procedure to evaluate nested multiple integrals with 
integrands being either a cosine series or a sine series 

Suppose that, up to an order A which depends on N± and 
N 2 , all the normalized QDD7V!,Af 2 coefficients a\ n ...\ 1 [multi- 
ple nested integral Eq. (01)] can be reduced to a single integral 
as either 



E / sin[ P6]d6 or 
E / cos[ P6}d6 



(48) 
(49) 



where we omit prefactors for simplicity. We shall show in the 
following subsections that this form arises in the evaluation of 
the QDDa^a^ coefficients. 

Moreover, we shall show in the following subsections that 
all a\ n ...\ 1= z coefficients with order n < N are of the form 
of Eq. (09), and hence that ax„--Ai=z vanishes since 



E sin[P0]|S = O 



(50) 



after performing the last integral. Therefore, the dephasing 
errors oz can be eliminated at least up to a remaining error of 

0(T N+1 ). 

We first note that regardless of the integration limits, all 
ax„ - Ai coefficients can be viewed as one integral nested with 
one order lower (n — 1 th order) coefficient, 



a A „... Al = / dOnhMall^ (51) 



where the superscript 0„ indicates that the upper integration 
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limitofa^" ^ is6> n ,not7r. Now assume that all the n—l th 
order coefficients a\ n _ 1 ...\ 1 are of the form of Eq. (|48l ) or 
Eq. (f49b . Then one could just proceed to the next order by 
substituting Eq. (l48t or Eq. (|49| i (with upper integration limits 
7r replaced by 9 n for the n — 1 th order coefficient o,\ _ 1 ...x 1 ) 
into Eq. (IBTl i. with the n-level nested integral having been re- 
duced to a two-fold nested integral. Therefore, under the as- 
sumption above, two-fold nested integrals are the basic units 
for evaluating multiple nested integrals. 

It follows from the result in the previous subsection that 
all /a„ are sums of purely sine or purely cosine functions. 
Combining this with Eq. d5TT l and the assumption that all the 
n — l th order coefficients a\„_i "Ai are of the form of Eq. ( |48T > 
or Eq. d49i i. there are only four possible types of two-fold 
nested integrals, which are presented on the left hand sides 
of Eq. d52l l. The results, on the right, follow simply from eval- 
uation of the 9 n -i integrals, followed by application of the 
product-to-sum trigonometric formulas ( 138a| >-( [38cl ). 



d9 n sm[p s 9 n ] c$ n _i cos[P c 0„_x] 
Jo 

J d9 n cos[( Ps ± P c )6 n ] 

d9 n cos[p c 9 n ] / d6»„_i sin[P s 6>„_i] 
Jo 

J dQ n cos[(p c ± P s )0 n ] 
d9 n cos[p c 9„] / <26> n _i eos[P c 0„_i] 



dO n sin[(p c ±P C )9„] 



(52a) 



(52b) 



(52c) 



J d9 n sm[p s 9 n ] J d#„-i sin[P,0„_i] ~ 

d9 n sm[( Ps ±P s )9 n } (52d) 



where the ± symbol is shorthand for, e.g., J d9 n cos[( Ps ± 
P c )9 n ] = / dO n cos[(p s + P c )9 n ] + J d9 n cos[( Ps - P c )9 n ], 
and where we have omitted irrelevant prefactors in front of all 
integrals. 

Note that the cos integrands on the right hand side of 
Eq. d52b . will yield 1 if their arguments happen to vanish. This 
conflicts with the requirement of Eq. d49l ). and would prevent 
us from proving that a\ n ...\ 1= z vanishes. Likewise, in order 
to proceed to the next order, say order n, none of the n — 1 th 
order coefficients a e x n _ . , in Eq. ( T5TT > may contain constant 
terms when expressed as a single integral of a cosine series. 
The reason that a constant is problematic is that it behaves 
differently from a cosine function under integration. The in- 
tegral of a cosine function with non-zero argument gives rise 
to a sine function, but the integral of a constant gives rise to a 
linear function. Therefore, if the integrand is a cosine series 
including a constant term, then after integration the result will 
not be a pure sine series any more. Furthermore, the problem 
cannot be resolved by carrying out the next integral. On the 



other hand, one need not worry about sine functions because 
sine functions with arbitrary angles will always result in co- 
sine functions after integration. 

Therefore, proceeding from n — 1 th order to n th order, sup- 
pose none of n — 1 th order coefficients a^" A contain 
a constant term. From Eq. ( l52b , due to the product-to-sum 
trigonometric formula, the problematic constant term will be 
generated when the new resulting argument p s ± P c in the 
cosine functions happens to vanish. When this happens to 
any one of the n th order coefficients aA„~Ai> there is no ad- 
vantage, when using our proof method, in proceeding to the 
n + 1 th order; this order is where the cosine arguments may 
start to be zero, and hence it sets a lower bound on the sup- 
pression order of the pure dephasing term. 



C. The suppression ability of the outer X -type UDD n 2 
sequence when JVi is even 

Let us define four function types we shall encounter in our 
proof. 



Definition 2. c" dd , c" ven , Q ven , and Q dd function types. Let 
k, q G Z with k arbitrary and \q\ < n. 

A c™ dd -type function is an arbitrary linear combination of 
cos[(2A: + 1)(JV 2 + 1)6< + qO] terms. 

A Cg Ven -type function is an arbitrary linear combination of 
cos[2k(N 2 + 1)8 + qff\ terms. 

A Cevcn~tyP e function is an arbitrary linear combination of 
sin[2fc(7V 2 + 1)6* + qO] terms. 

A C^ dd -type function is an arbitrary linear combination of 
sin[(2fc + 1)(JV 2 + 1)6* + qO] terms. 

Whenn < N 2 we have (2k+l)(N 2 +l)+q ^ 0. Therefore, 
by definition, all c" dd -type functions will in this case have no 
constant 1 (the problematic term). The c™ vcn -type functions 
are allowed to have a constant term. 

From Eqs. (I42l-d45ll. for even inner decoupling order N±, 
there are only two kinds of integrands: fo and f x are Civen" 
type functions while f z and f y are c^ dd -type functions which, 
as we just remarked, do not have the constant 1 term. There- 
fore, it immediately follows from Eqs. (|49| i and (l50l l that the 
first order normalized QDDa^ jv 2 coefficients az = J Q n f z d9 
and ay = f~ f y d9 vanish. 

Next, let us consider the second order terms (two-fold 
nested integrals), as in Eq. (l52l . We introduce a binary oper- 
ation which (1) evaluates the first integrand, (2) multiplies 
the outcome with the second integrand, (3) applies the appro- 
priate product-to-sum trigonometric formula. Substituting the 



Ddd 



or Civen"tyP e functions into Eq. ( 1521 . we then have 



odd 



c 1 

^even 

r 1 G)^ 
odd w Seven 



-odd 



-"odd 



odd © c odd 
©Ce, 



c 1 

Seven 



C 2 

Seven 
even Seven 



(53a) 
(53b) 
(53c) 
(53d) 



where we omitted the second integration symbol. 
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If we disregard the n superscript of c" dd and C" V cn m 
Eq. ( I5"3l ). the set {Coven, c dd} constitutes the abelian group 
Z2 under the binary operation 0, with the identity element 

Seven- 

On the other hand, the superscript of the resulting function, 
c odd or Ceven m Eq. d53l ), is just the sum of the superscripts of 
the first and second integrands (c^ dd or Coven)- Accordingly, 
the binary operation acts as integer addition for the super- 
script 11. 

Let us now consider the ?i-fold nested integral implied by 
Eq. ( T5Tb . Because of the closure property of the group Z%, 
integer addition of the superscripts n of c Q l dd , and CTven' an d 
the fact that no c™ dd -type function with n < N2 contains the 
constant 1, we can conclude that such an n-fold nested integral 
with n < N2 and with each integrand being either c^-type 
or Ce Ven -type functions can be reduced to be either J c o l dd d0 n 
orjQ dd de n . 

More specifically, note that in Eq. d53l > the first two lines 
have an odd number of c^ dd functions and result in c^ dd , while 
the last two lines have an even number of c* dd functions and 
result in Ccvon- When we continue the nesting process using 
these rules, the odd or even property is maintained while the n 
superscript grows by one unit each time. In other words, due 
to Z2 group multiplication rules [Eq. d53l l without the super- 
scripts n], we have the following lemma: 

Lemma 5. Provided n < N2, all n-fold nested integrals, with 
each integrand being either a c^ dd -type or C^ ven -type function, 
can be written as 

1. J c™ dd d9 n if there is an odd total number of c\ dd -type 
integrands in the n-fold nested integral, 

2. J Ceven if there is an even total number of c^ dd -type 
integrands in the n-fold nested integral. 

Next, let us determine the parity of the number of c^-type 
functions appearing in the QDD^TVa coefficients. Consider, 
e.g., a A „...Ai=z- Recall that f and J x are Ceven-tyP 6 func- 
tions while f z and f y are c odd -type functions. Consulting the 
last column of Part 1 of Table Hill (the second a\ n ...\ 1= z col- 
umn), we see that there is an odd number of f x (Civen) an d fy 
( c odd)' an d an even num ber of f z (cj dd ). Therefore, in this 
case, we have an odd+even=odd number of c^-type func- 
tions in a\ n ...\ 1= z- Similarly, consulting all other columns of 
Part 1 of Table [TlTl it turns out that all possible combinations 
generating a\ n ...\ 1= z or a\ n ...\ 1= Y contain an odd number 
of CQ dd -type functions. It now follows from Lemma [5] and 
then Eq. d50b that all a\ n ...\ 1= z = a\ n -\ 1= y = if the or- 
der n < N2. [Note that this counting argument is unaffected 
by the move from / to /, since this move was due to a change 
of integration variables — see Eq. (1411 .1 

In conclusion, the inner UDD^ sequences with even order 
Ni do not affect the suppression effect of the outer UDD jy 2 
sequence, i.e., the az error is always removed up to the ex- 
pected order N2 when the order JVi of the inner sequence is 
even. This proves the first row of the az part of Table IIVI 

In addition, we have just shown that when the inner order 
N\ is even, the outer X-type UDDjy 2 sequence also elimi- 
nates the ay error up to the outer decoupling order N%. Since 



we have shown in Sec. [Ill] that the ay error is suppressed to 
order N\ when N% is even, or Ni + 1 when N2 is odd, one 
can conclude that when the inner order Ni is even, ay is sup- 
pressed to order max[iVi , N2] when N2 is even, and to order 
max[7Vi + 1, N2] when N2 is odd. This completes the proof 
of the ay part in Table IIVI 



D. The suppression ability of the outer A'-type UDDjv 2 
sequence when Ni is odd 

The main difference between the analysis in this subsection 
and the previous one is that f x and f y are interchanged in 
terms of which function is cosine or sine — see Eqs. (I44b-(l47b. 

Also, note that, from Eqs. <@3, d43j, (g6]l, and ATI ), for odd 
inner decoupling order JVi, fo is a Ceven _t yP e function, f z is a 
Codd-type function, f x is a ci vcn -type, and /„ is a Co^d-^P 6 - 
We shall use these facts throughout this subsection. 

Our procedure is to start from the first order QDD coef- 
ficients, then the second order, and finally the general, nth 
order. 



1. The first order terms a\ 1 

It immediately follows from the fact that f z is a c^ dd -type 
function and from Eq. 69) that a z = £ f z d9 = 0. It also 
immediately follows from the function types that fo and f y 
are of the form of Eq. ( f48b . The only function that deserves 
special attention is f x . 

As discussed in Sec. IIV Bl in order to proceed to second 
order, none of the modulation functions can contain a constant 
1 term. However, Definition [2] allows cj, vcn -type functions to 
have such a term. Accordingly, before applying Eq. ( BTt to 
the second order case, we should check whether f x {0) has a 
constant term. 

Suppose f x (6) has a constant 1 term and then separate the 
constant 1 from the other cosine functions with non-zero ar- 
guments as follows, 

f x {6) =^d p cos[p#]+r (54) 

where p = 2k(N2 + 1) ± q with \q\ < 1 an integer, r a coef- 
ficient of the constant 1 term, and d p coefficients of cos[p6>]. 
Then the first order normalized QDD jvj n 2 coefficient of the 
ax error reads 

ax = f h{d)d9 
Jo 

= Y,dvMpO\\e=Z + rO\ B eZZ 

p=£0 

= + rir (55) 

Now, since we already proved in Section lHII that a,\ n \ 1= x — 
for n < Ni, and since in the first order case 77 = 1 and 
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hence n < Ni always holds, it follows that r = 0. Therefore, 
f x does not contain a constant 1 term. 

In summary, now that we have shown that a\ 1= z = and 
that all the first order normalized QDD^ jv 2 coefficients a\ 1 
are of the form of either Eq. d48t or Eq. d49l ), we can proceed 
to the second order case. 



2. The second order terms a\ 2 \ 1 



the integrand c^ vcn does not have a constant 1 term if N\ > 2. 

In order to proceed to the next order none of the integrands 
may contain a constant. Therefore, our results show that if 
Ni,N2 > 2, one can indeed proceed to the next order. On 
the other hand, if Ni = N2 = 1 we can only conclude that 
a\ 1= z = 0, while if N% = 1 and N2 > 2, we can only 
conclude that a\ 1= z = a\ 2 \ 1 =z = 0, but not that the third or 
higher order Z-type QDD coefficients are zero. 



From Eq. dBll l. the second order normalized QDD n 1 jv 2 co- 
efficients is 



3. The n 1 order terms a\ n 



d6 2 fx 2 (6 2 )al 



(56) 



After additionally applying the trigonometric product-to-sum 
transformation, and using the operation defined above 
Eq. (l53l , we can write 



OAaAi = OA 2 a A 2 



(57) 



Next, in Eq. d57K let us substitute Ceven mt0 tne integrand 
of ai, cl vcn into a x , Codd int0 a Y, and c odd into az- The 
resulting set of all <ja 2 Ai can be arranged into a multiplication 
table, Table I VII where the entries in the top row are the types 
of the first integrand and the entries in the left-most column 
are the types of the second integrand. The remaining entries 
are the results of applying the binary operation between the 
elements of the first row and column. 

From Table IVI1 the superscript of the resulting func- 
tion is again the sum of the superscripts of the first and 
second integrands. Hence the binary operation again 
acts as integer addition for the superscripts. Moreover, 
disregarding the superscripts n, Table [VI] shows that the 
set {Cevcn, c C vcn, Codd, c dd} forms the abelian Klein four- 
group, i.e., the Z2 x Z 2 group, under the binary operation 0. 
The key observation from Table |Vl]is that the algebra of the 
subscripts A2A1 of aa 2 Ai works as the Pauli algebra without 
the anti-commutativity property, which is isomorphic to the 
Klein four-group algebra by mapping the identity / to Ccvcn, 
X to c cvcn , Y to (odd, and Z to c od d- 

Accordingly, the results of Table |Vl]can be summarized as 
follows, 



a A 2 Ai = Z 

a\ 2 \ 1= i 



Codd d0 



d9 



Codd d0 



So 



de. 



(58a) 
(58b) 
(58c) 
(58d) 



We can conclude that a\ 2 \ 1= z — if > 2, since then (by 
definition) c^ dd does not contain a constant 1 term. 

We have already proved in Section|III]that a\ 2 \ 1= x = if 
Ai > 2. By the same argument as Eq. (1331 ), this implies that 



The procedure we described for the first and second orders 
applies to higher orders, until one reaches the order A where 
some resulting integrands begin to include constant 1 terms. 

To obtain the nth order QDDjVj^ coefficients we proceed 
by induction on n. We have already established the case of 
n = 1 and n — 2. Let us assume that 



aA„---Ai = Z 

a\ n ...\ 1 =x 
a\ n ...\ 1= Y 
a\ n ...\ 1= i 



Cid^ 



Q dd de 



r 71 

/ Coven dO. 
JO 



(59a) 
(59b) 
(59c) 
(59d) 



where none of these integrals contains a constant 1 term in 
their integrand, and prove that the same integrand form holds 
for n + 1 (but not necessarily that there is no constant 1). 
Indeed, using the definition of the operation and Eq. ( BTT ). 
we have 



a A, l+1 ©a£... Ai 



(60) 



Due to the induction assumption [Eq. (|59l l the situation is 
now identical to the one we analyzed for n = 2, in particular 
inEq. d58l ). Therefore Eq. d59l holds with n replaced by n+1. 

This can also be understood without induction as being due 
to the isomorphism between the set {Ceven, c eV en, Codd, c dd} 
and the set {/, X, Y, Z} (the subscripts of a\ n ...\ 1 ), and the 
addition of superscripts under the operation. 

To figure out up to which order A Eq. ( 1591 holds, one must 
examine when the c™ dd or c^^-type functions begin to have 
constant 1 terms. The c" dd -type functions will by definition 
not contain constant 1 terms until order n = N2 + 1. On the 
other hand, due to Eq. ( 159bl i and a\ n ■■■\ 1 =x = for n < N\ 
(proven in Sec. [Hi]), c™ V en m Eq. ( I59bb is guaranteed to have no 
constant 1 term until order n = N\ + 1, by a similar argument 
as that leading to Eq. ( 1331 ). In conclusion, 



Lemma 6. For QDD^^ with odd N\, all nth order nor- 
malized QDDn 1 n 2 coefficients a\ n ...\ 1 with n < min[Ai + 
1, A 2 + 1] can be written as Eq. J591 , and none of the in- 
tegrands in Eq. d59t contain a constant 1 term when n < 
min[Ai, A 2 1. 
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a / • Ccvcn 


a X • C C vcn 


a Y ■ Codd 


a z ■ c odd 


a I ■ Ccvcn 


ail=l ■ Ccvcn 


a IX=X ■ c cvcn 


aiY=Y ■ Codd 


aiz=z ■ c z odd 


"*X ■ ^cvcn 


a XI=X ■ c cvcn 


aXX=I ■ Ccvcn 


axY=z ■ c z odd 


axz=Y ■ Codd 


a Y ■ Codd 


ay/=F : Codd 


ayx=z ■■ c 2 odd 


a Y Y=I ■ Ccvcn 


a Y Z=X ■ Co Von 


az ■ c odd 


azi=z ■ c^ dd 


azx=Y ■ Codd 


azY=x ■ c^ vcn 


0-ZZ=I ■ Ccvcn 



TABLE VI: The group structure associated with the second order QDD^ 1 n 2 coefficient a\ 2 \ 1 = a\ 2 © . 



It follows immediately from Lemma|6]and Eq. ( l50b that the 
first min[JVi , iVa] orders of a\ n ...\ 1= z vanish. However, in 
fact we can show more, namely that a\ n ... \ 1= z = for all 
n < min[7Vi + 1, N%]. Suppose that ATi < 7V2 and consider 
the special case n = Ni + 1. In this case it follows from 
Lemma[5]thataA Wl+1 ...Ai=x = Jo* c odd~ 1 ^> tne argument of 
the function c^ 1 is (2fc+l)(A^ 2 + l)^+g6>, and \q\ < iVi+1. 
Since iVi < -/V2 this argument cannot vanish, and it follows 
thataA N \i=z = 0. In conclusion, ax„— Ai=z = for all 
n < min[A r i + 1, N2}, which proves the last row in Table IIVI 

In summary, if the inner decoupling order N\ is odd and 
N2 < Ni + 1, the outer UDDjv 2 sequence always suppresses 
the dephasing error Z to the expected decoupling order N2, as 
then minfA^i + 1, N%] = N2. In contrast, if the inner decou- 
pling order N\ is odd and N2 > Ni + 1, the outer UDDjy 2 
sequence suppresses the dephasing error Z (at least) up to or- 
der Ni + 1, which may be smaller than the expected outer 
decoupling order N2 . Thus, if the order of inner level UDD n 1 
sequence is odd, this may inhibit the suppression ability of the 
outer UDD jv 2 sequence. 



V. COMPARISON BETWEEN OUR THEORETICAL 
BOUNDS AND NUMERICAL RESULTS 

In Ref. Ir42tl the QDD sequence was analyzed numerically 
and the scaling of the single-axis errors was determined on the 
basis of simulations, for Ni and N2 in the range {1, . . . , 24}. 
These simulations are in complete agreement with our ana- 
lytically bounds for n x and n y , as given in Table [IV] They 
are also in complete agreement with our bound for n z when 
N\ is even. Thus we can conclude that it is likely that our 
bounds are in fact tight in these cases. There is, however, 
one discrepancy: when Ni is odd our analytical bound yields 
n z = min[iVi + ljiVa], while the numerical result found in 
Ref. Ir42tl is n z = min[27Vi + 1, N2]. Thus, in this case our 
bound is not tight. We attribute this to the fact that the method 
we used in Sec. [IV] does not use the full information con- 
tained in the integrands, i.e., we discard all Fourier coeffi- 
cients. Specifically, if a\ n ,,,\ 1 —z contains a constant term, 
namely, cos [P6] with P = 0, or does not end up in the form 
of Eq. d49l i. it is still possible that a\ n ...\ 1= z vanishes because 
a sum of non-zero terms could be zero when combined with 
the right Fourier coefficients. Thus, our method of analysis 
merely yields a lower bound on the decoupling order of the 
pure dephasing error. It is an interesting open problem to try 
to improve this bound so that it matches the numerical results 
of Ref. JH. 



VI. SUMMARY AND CONCLUSIONS 

The QDD sequence, introduced in Ref. fiUl . is, to date, 
the most efficient pulse sequence known for suppression of 
single-qubit decoherence. In this work we provided a com- 
plete proof of the validity of this sequence, i.e., we proved its 
universality (independence of details of the environment) and 
performance. Our work complements an earlier proof Ir44ll . 
which was restricted to even order inner UDD sequences. 
However, our results go beyond a validity proof of QDD. For, 
in this work we also elucidated the dependence of single-axis 
error suppression on the orders Ni and N2 of the inner X-type 
and outer Z-type UDD sequences comprising QDOat^, re- 
spectively. Our results are stated in Theorem 1 . Let us briefly 
summarize our method and main findings. 

Our general proof idea was to analyze the conditions un- 
der which, for each error type a\, the nth order ODD^r^ 
coefficients [Eq. ( TT8b l vanish. We used two complementary 
methods. In the first method, we expressed the QDD coeffi- 
cients ax n ...Ai m terms of UDD coefficients by splitting each 
of OA n —Ai' s nested integrals into a sum of sub-integrals over 
normalized outer intervals. We were then able to conclude 
that a\ n ...\ 1= x and <XA n ...Ai=Y vanish when n < Ni due to 
the vanishing of the UDDa^ contributions. For the ay error, 
still as part of the first method, we showed that an additional 
order vanishes due to a parity cancellation effect involving the 
outer sequence. However, this additional cancellation cannot 
be attributed to the vanishing of a corresponding UDD co- 
efficient. In the second method we considered the case of 
a\ n — Xi=z, for which we provided an analysis based on the 
evaluation of integrals of trigonometric functions. We showed 
that their properties under nested integration can be mapped to 
the Abelian groups Z2 (for even Ni) and Z2 x Z2 (for odd N\). 
Using this we provided a proof by induction for the vanishing 
of o-x n — \!=z> and, when Ni is even, also for aA„...Ai=r- 

The overall summary of our results is that aA„~Ai=A = 
Vrt < N, where N is the decoupling order given in the last 
column of Table IIVI We now provide a recap of these results, 
including a semi-intuitive explanation based on the idea of in- 
terference between the modulation functions. 

Starting from the simplest case, we showed explicitly that 
independently of the order of the outer X-type sequence, the 
inner Z-type UDDjVi sequence always achieves its expected 
error suppression order, i.e., the <tx and ay errors are sup- 
pressed to the inner decoupling order N±. Since ax errors 
commute with the pulses of the outer sequence they are not 
suppressed any further. 

The story is more complicated for the ay and az errors, as 
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they are both suppressed by the outer sequence. 

For the cry error, the parities of the inner and outer sequence 
orders cause the decoupling order to vary between Ni, Ni + 1, 
and N2- Consider first the even N2 case. An intuitive expla- 
nation for the corresponding parity effects is the following. 
For even N\ , the modulation functions f y and f z are in phase, 
namely both have a sin[(2/j + 1)(7V2 + 1)6*] dependence [recall 
Eqs. ( f33T > and (|39]ll. The outer X-type UDDtv 2 sequence, with 
its f z modulation function, is then fully effective at eliminat- 
ing the cry error, with the result that oy is eliminated to the 
expected decoupling order max[Xi, X2]. However, when N\ 
is odd, f y has a cos[(2fc + 1)(N 2 + 1)0] dependence, which 
is 90 degrees out of phase with f z . In this case f y and f z in- 
terfere destructively with one another, and the outer sequence 
does not help to further suppress ay. The result is that cry is 
only eliminated to order N\ . 

Now consider the case of odd N2. This case gives rise to 
the anomalous N% + 1 suppression order. The reason is that 
when N2 is odd, the modulation function f z is odd with re- 
spect to the midpoint of the total sequence duration, while f x 
and f y are both even. It is this oddness of the outer sequence 
modulation function (f z ) which helps to suppress the error cry 
to one more order, due to a cancellation of terms with equal 
magnitude but opposite sign [Eq. d32ll. This gives rise to a 
cancellation to order max[Xi + 1, N2] when Ni is even and 
the inner sequence does not interfere with the outer sequence, 
or to order Ni + 1 when Ni is odd and the inner sequence 
does interfere with the outer sequence. 

Thus, suppose we fix N2 so that it is even (odd) and greater 
than Ni (Ni + 1). We should then see the suppression order 
of <Ty switch between Ni (Ni + 1) and N2, as Ni is increased 
from 1 to N2, a phenomenon which was indeed observed in 
the numerical simulations of Ref. l42tl . 

If the inner order Ni is even, the outer X-type UDDat 2 se- 
quence always suppresses <rz to the expected decoupling or- 
der N2- This has the same intuitive origin as the cry case. 
Namely, for even Ni, f y and f z are in phase, i.e., both have a 
sin[(2fc+l)(A r 2 + l)0] dependence, and so are able to suppress 
o~z to the expected order. However, when Ni is odd, the de- 
pendence of f y is cos[(2fc+ l)(N2 + 1)9, which is 90 degrees 
out of phase with f z . Therefore again f y and f z interference 
destructively, and the outer sequence does not suppress the 
error oz to the expected order. 

In more detail, if the inner order Ni is odd and N2 > 
Ni + 1, our proof method shows that the outer X-type UDD jv 2 
sequence suppresses the oz error at least to order Ni + 1, 
which is less than the expected outer decoupling order X2. 
Hence, if this lower bound is saturated, one can see a satu- 
ration effect in the decoupling order of oz, which starts at 
N2 = Ni + 2 when we fix odd Ni and increase N%. Thus, odd 
Ni can hinder the suppression ability of the outer sequence. 

The numerical results of Ref. f42ll confirm that odd Ni can 
hinder the suppression ability of the outer X-type UDDjv 2 se- 
quence. However, the actual saturation effect in the decou- 
pling order of uz begins at N2 = 2Ni + 2, higher than our 
lower bound of N2 = Ni + 2. A new method may be needed 
to explain the remaining vanishing orders from A^i + 2 to 
2Ni + 1. 



The inhibitory effect of odd inner decoupling order Ni dis- 
appears when Ni is large enough. Specifically, whenever 
Ni > N2 — 1 the outer X-type UDD^ sequence suppresses 
(Tz to the expected decoupling order N%. This makes intu- 
itive sense because when Ni is large enough the outer X-type 
UDD at 2 sequence "views" the effective Hamiltonian resulting 
from the inner Z-type UDD^ sequence — which has time de- 
pendence 0(T Nl+1 ~ N2 ) — as time-independent relative to its 
"error cancellation power" 0(T N2 ). 

Despite this complicated interplay between Ni and N2, our 
proof yields the simple result that the QDDjVjAT;, sequence 
suppresses all single-qubit errors to an order > min[Xi, N2). 
This matches the numerical results in l42tl . so that our bounds 
appear to be optimal in this regard. We conclude that to attain 
the highest order decoupling from the QDDn^^^ sequence 
with ideal, zero-width pulses, one should use either an even 
order inner UDD sequence, or ensure that Ni > N2 — 1 if Ni 
is odd. 

A natural generalization of the work presented here is to 
NUDD with different sequence orders IHfll . We look forward 
to experimental tests of the properties of the QDDtvj^ pulse 
sequence predicted in this work. 

Note added: After this work was completed and while it 
was being written up for publication we became aware of 
a different, elegant proof of the universality of NUDD and 
in particular QDD 14711 . Our approach differs not only in 
methodology but also in providing a complete analysis of the 
single-axis errors. 
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Appendix A: The form of a\ n . after outer interval 
decomposition 

We shall derive Eq. d22b by splitting each integral of 
a\ n ...\ l [Eq. (fSTJ] into a sum of sub-integrals over the normal- 
ized outer intervals Sj in Eq. (f3ab . Since a\ rl ...\ 1 comprises a 
series of time-ordered, nested integrals, our procedure for de- 
composing a\ „...Ai is to split its nested integrals one by one, 
from to T)' 1 '. 

We call the sub-integral over the jth outer interval "sub- 
integral-j". Suppose the integral of the integration variable 
r/W follows the sub-integral-j^ +1 ' of the previous variable 
rj^ +1 \ By splitting the integral of 7/^ with respect to the 
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normalized outer intervals and using Eq. ( fTST ). we have 



Appendix B: Time symmetry of the UDD pulse intervals 



e /ftC? w ) p w (Aia) 



+f Pl (j (e+1) ) / yWW 5 (Mb) 
.1 J (f+1) -1 

/cM^Jd^/ftU^)*^!) (Aid) 

To obtain Eqs. ( lAlcb and ( lAldb we rescaled fa t (ff^) in 
Eq. (lAlab and ( lAlbt individually with 



(A2) 



for each outer interval s j(f ) , thus obtaining f at (fj^'). In this 
manner fativ ) is me same function for all the outer inter- 
vals, so that L f ae {rf 1 ' ) dff^ can be taken out from the sum- 
mation, as shown in Eqs. ( I A let and ( lAldb . 

Recall the time-ordering condition, rf n > > rf 71 ^ 1 ' > 
■ ■ ■ V > ^ ■ It has a consequence that in Eq. dAldl ), sub- 
integrals over any two adjacent variables r/^ and rj^ +1 ^ are 
nested, as they are in the same outer interval, number 
In this case jyW < rf +1 > . 

In contrast, if the sub-integrals are in different outer inter- 
vals (automatically time-ordered), then the sub-integral over 
is not nested inside the subintegral over rf i+1 >, but inte- 
grated over its entire outer interval independently, as in Eq. 
( Kiel . 

Let rt = * denote the time-ordering of outer intervals as 
in Eq. (lAlcl ). and let 7^ = denote the integral time-ordering 
inside a given outer interval as in Eq. dAldl i. Accordingly, 77 
describes the relation between the adjacent variables 
and tjW. 

As we have just shown, each integral of (Za„...Ai can always 
be split into two parts, Eq. ( lAlcb and ( lAldb . with one excep- 
tion: if j" +1 > = 1, the subsequent sub-integral of variables 
rf® will only contain the term Eq. dAldl ). Moreover, both 
Eq. dAlcl ) and dAldb contain an effective inner part (the part 
that depends on f at ) and an effective outer part (the part that 
depends on fp t ). Therefore, by substituting Eqs. dAlcl ) and 
( lAldb into each integral of aA n ...Ai, m sequence from rj( n > to 



(i) 



gsa„...Ai can ^> e written as an inner part $ ln [Eq. d23l ll 
multiplying an outer part $ out [Eq. (f25])] over all the pos- 
sible integration and summation configurations. Each such 
configuration can be denoted by an ordered set of symbols 
{r„_ 1 r n _ 2 . . • ri}. Thereby, we obtain Eq. (1221 as the repre- 
sentation of a\ n ...Ai after this decomposition. 



Due to the identity sin 9 
Eq. d3ab satisfies 



sm 7T 



r (2j-l)7r i 
i 2(JV+1) 1 



in s, 



sm 



(2j - 1)tt 
2(JV+1) 



sin 



(2j - 1)tt 
2(N + 1) 1 
(2N + 2 - 2j + 1)tt 
2(N + 1) J 
= sin[ (2(^ + 2-j)-l) 7 r 
1 2(7V + 1) J 

Therefore, we have proved that Sj = sjy+2-j [Eq- OTbl, 
which shows that the UDD pulse intervals are time symmet- 
ric. There is, however, a difference between even and odd N: 
when N is odd every interval to the left of center is paired 
with an interval to the right of center. When N is even the 
central interval is unpaired. E.g., for N = 1 we have two, 
paired intervals: si = S2. When TV = 2 we have two paired 
intervals, si = S3, and an unpaired interval S2. 



Appendix C: The parity of the inner order Ni determines the 
parity of f x 

Since the inner pulse sequences under the piecewise linear 
variable transformation Eq. (f33j) still have the UDD jvj struc- 
ture, the rescaled inner pulse intervals remain time symmetric: 



7j,N x +2-k — ty,iVi+l-fc- 



(CI) 



When the inner decoupling order Ni is even, the parities of 
Ni + 2 — k and k are the same, so that 

f (a\ _ H -1 )*" 1 6 e [9j,k-l>Qj,k) fro . 

Hence f x (9) is even inside each outer interval. 
When the inner decoupling order N\ is odd, 



fx(8) 



(- 1 ) fe " 1 2 k (C3) 



where the sign difference between the second lines of Eq. ( IC2b 
and Eq. dC3l ) arises from the opposite parities of JVi + 2 — fc 
and fc. Accordingly, f x (9) is odd inside each outer interval. 

Note that the sequence of rescaled inner intervals 
{^•.fclfii 4 is repeated for all values of j e {1, . . . , N 2 + 1}. 
As a result the three modulation functions f x (6), f y (9), f z {9) 
are periodic, with respective periods N ^ +1 , jv^+i ' JV2+1 ' ^ n 
this sense, the variable transformation r\ = sin 2 (9/2) intro- 



duced in 13011 . which emerges naturally from the time structure 
of UDD sequence Eq. ©, is unsuitable for our QDD proof. 
The reason is that despite the fact that the outer X-type pulses 
intervals are rescaled to be equal, the timing patterns of the 
inner sequences in different outer intervals are no longer the 
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same. 



Appendix D: Fourier expansions of G(6) 



G(9) in Eq. J40bt takes the following form up to a multi- 
plicative constant: G{6) = Sj, where 9 £ [ , j^tx)- The 
symmetry property (f3jj) implies that G(6) can be written as 



G(0) = X>siiU 



(Dl) 



£=1 

Let us now compute the expansion coefficients: 

91 = [ G{6)sml6d6 
n/ 2 Jo 

N+i r e< 



7T 



Sj I sin IBdB 



j = l -">j-l 



N+l 



= ; Sin ■ 



— Y 



sm 



(2j - 1)tt 

tt^ 2(iV + 1) 2(AT + 1) 

(cos — cos £0j _ i ) 

' • (2j-l)7r 



AT+1 



sm — -f '— x 

tt^ 2(7V+1) ^ 2(N + 1) 



(-2)sin£ 



smf 



(D2) 



2 2 

where we used the sum-to product formula in the third equal- 
ity. Due to 2 3 ~ 1 = 2 {n+i) an< ^ me product-to sum for- 
mula, we have 

4 7r £ir 

gi = -£ s[n WTT) s ' m WTT) x 

N+l 



. (2j - 1)tt . (2j-l)7r 
> sm — sm £ — 

^ 2(/V + l) 2(JV+1) 

4 . 7T . t-K 

— sm — sm — r x 

Tit 2{N + 1) 2{N+\) 

1 iV+1 

— cos t^—^ — cos ■ 

2 ^ 2(N + 1 
j=i 



(D3) 



(*-l)(2j-l)7r _ (£ + 1)(2j-1)tt 

2(JV+1) 



Considering the sum over j we have 



2(N + 1) 



j'=i 

^±1)J7T (£±1)7T 



E 



cos 



(iV+1) 2(7V + 1) J 



N+l 



■ ( t I 1 j TT ^ ^ - I -C ZC J- J J 77 

Re[e" l2 ( JV + 1 > \ e l ] 
i=i 

■ (g±i)x (f± i )jr j _ e »(i±i)fl- 
Re[e ^(iv+Dg' N+i (f±i)7r ] 

1 - e l 

l-cos(^±l)7r-isin(£±l)7r 

■ (l±l)ir" ■ (t±l)7T J 

g S 2(JV + 1) _ g*2(]V + l) 

sin(i ± 1)tt 



2 sin 



2(W+1) 



(D4) 



where in the third equality we used the geometric series for- 
mula. The last expression vanishes if I ^ 2k(N + 1) =p 1. The 
only values of I for which gg does not vanish are 2k(N+l)^fl. 
Therefore, finally 



G{6) = E gk >i S H 2 KN + 1)0 + qO]. (D5) 

fe=0 q=±l 
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